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Abstract 

Topical phenomena in high-energy physics related to collision experiments of heavy nuclei ("Lit- 
tle Bang") and early universe cosmology ("Big Bang") involve far-from-equilibrium dynamics 
described by quantum field theory. One example concerns the role of plasma instabilities for 
the process of thermalization in heavy-ion collisions. The reheating of the early universe after 
inflation may exhibit rather similar phenomena following a tachyonic or parametric resonance 
instability. Certain universal aspects associated to nonthermal fixed points even quantitatively 
agree, and considering these phenomena from a common perspective can be fruitful. 



1. Introduction 

Collision experiments of heavy nuclei at the Relativistic Heavy Ion Collider (BNL), 
the Large Hadron Collider (CERN) and future experiments at the FAIR facility (GSI) 
involve far-from-equilibrium dynamics for strongly interacting matter described by quan- 
tum chromodynamics (QCD). Available data from RHIC reveals remarkable, unexpected 
properties such as rapid apparent thermalization with robust collective phenomena. 
There seems to be a good description of a wide range of data using hydrodynamical 
models with near-perfect fluidity, provided one assumes that suitable initial conditions 
form in about 0.5-2 fm/c |l|2j . The explanation of these findings from QCD provides a 
challenge for theory. Strong interactions can play an important role for QCD at acces- 
sible energies. However, there are theoretical indications that essential nonequilibrium 
dynamics could affect the picture dramatically. Complete local thermal equilibrium may 
not be necessary for the application of hydrodynamics |3|4j . For sufficiently high energy 
densities nonequilibrium instabilities in anisotropic plasmas have been identified as the 
parametrically fastest processes, which may help to explain a rapid isotropization of the 
equation of state relevant for near-perfect fluid descriptions [514161718] . 

Nonequilibrium instabilities may also play an important role for the process of thermal- 
ization in the early universe after inflation |9ll0lllll2ll3ll4j . In inflationary cosmology, 
the universe at early times expands quasi-exponentially in a vacuum-like state. During 
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this stage of inflation, all energy is contained in a slowly evolving inflaton field. Eventu- 
ally the inflaton field decays and transfers all of its energy to particles, thereby starting 
the thermal history of the hot Friedmann cosmology. In chaotic inflationary models for 
a wide range of couplings the particle production from a coherently oscillating inflaton 
occurs in the nonperturbative regime of a parametric-resonance instability. This picture, 
with variation in its details, extends to other inflationary models. 

Though heavy-ion collisions and inflationary cosmology involve energy scales many 
orders of magnitude apart, their description requires similar quantum field theoretical 
techniques and certain universal aspects can even quantitatively agree. Instabilities lead 
to exponential growth of occupation numbers in long wavelength modes on time scales 
much shorter than the asymptotic thermal equilibration time. Though the underlying 
mechanisms are very different for QCD and for scalar inflaton dynamics, the subsequent 
evolution after an instability follows similar patterns: After a fast initial period of ex- 
ponential growth the dynamics slows down considerably. The system is still far from 
equilibrium at this stage and the subsequent slow evolution is characterized by power- 
law distributions. A new class of infrared scaling solutions, where a characteristic large 
correlation length leads to independence of long-distance properties from details of the 
underlying microscopic theory, has recently been found in the context of early-universe 
reheating dynamics |15| . The signatures can be dramatic. Nonthermal scaling solutions 
can exhibit strongly enhanced low- momentum fluctuations p^^, as compared to a ther- 
mal high-temperature distribution p^^. This long-distance behavior is distinct from 
scaling associated to (weak) turbulence, such as ~ which is relevant for character- 

istic properties at shorter distances [16 . It is remarkable that the same high-momentum 
scaling exponents known from inflaton dynamics have recently also been found to describe 
gluonic properties within lattice QCD in the classical-statistical gauge theory limit [17]. 
Research concerning quantitative agreements of far-from-equilibrium properties in heavy- 
ion collision experiments and early universe dynamics is still in its infancies. This article 
highlights some known aspects and important open questions. 



2. Collisions of heavy nuclei and QCD plasma instabilities 

A wide range of experimental RHIC data, including elliptic flow measurements for 
momenta of less than about 1-2 GeV, can be described using hydrodynamical models with 
near-perfect fluidity [112] . The essential assumption of ideal hydrodynamic descriptions 
is the presence of an almost constant equation of state P — P{e), which relates energy 
density e to pressure P given by a nearly diagonal stress tensor 

c PS., (1) 

in the local fluid rest frame. This is needed to close the system of equations obtained 
from the conservation of the energy-momentum tensor. Relation ([T]) is true for isotropic 
systems if typical relativistic excitations have random directions, even if their energy 
distribution is far from thermal. Let "T" denote the characteristic momentum of typical 
excitations out of equilibrium, which may be identified with the saturation scale Qs at 
time Qj^ in the saturation scenario 18J. The relevant dephasing time which is required 
to obtain an almost constant or "prethermalized" equation of state for isotropic systems 
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is tpt ^ 0{1/T) 3 . For sufficiently weafcly coupled theories this time scale is much faster 
than the inverse relaxation rate of arbitrarily small departures from equilibrium given by 



lax 



o(i//r) , 



(2) 



which is associated to large-angle scattering or near-collinear splitting processes among 
quarks and gluons in QCD with gauge coupling g [f9 . For anisotropic systems this 
relaxation time would characterize the time for isotropization in the absence of nonequi- 
librium instabilities [50]. However, in QCD nonabelian Weibel instabilities can drive local 
isotropization on a time scale, which is faster than ordinary perturbative scattering pro- 
cesses 4 . For sufficiently weak coupling and order-one anisotropy the inverse instability 
growth rate 



is parametrically faster than the inverse relaxation rate ^ [5]- Consequently, under- 
standing when a hydrodynamic description can first provide a good approximation is the 
same question as understanding the fastest characteristic process driving isotropization. 

Most quantitative estimates of these processes rely on a separation of scales between 
suitably defined 'soft' and 'hard' momenta for sufficiently small characteristic running 
gauge couphng [S]. In this case the gluon degrees of freedom are described as classical 
fields for soft and classical particles for hard scales. The respective hard- loop effective 
theory of soft excitations is a nonabelian version of the linearized Vlasov equations of 
traditional plasma physics, which are based on coUisionless kinetic theory for hard par- 
ticles coupled to a soft classical field[l] This approach neglects quantum corrections and 
may also be considered as an approximation of the classical-statistical field theory limit 
of the respective quantum gauge theory [7l8j . Classical-statistical lattice gauge theory 
provides a quantitative description in the presence of sufficiently large energy density or 
occupation numbers per mode. The simulations are done by numerical integration of the 
classical lattice equations of motion and Monte Carlo sampling of initial conditions. In 
the following we concentrate on classical-statistical lattice gauge theory simulations of 
Ref. ^ . The different approximations agree to a large extent in their predictions for char- 
acteristic growth rates at early times. They apparently disagree, in particular, concerning 
subsequent evolutions to turbulent flow which is described below. 

Collisions of heavy nuclei leave behind a plasma of quarks and mostly gluons in a 
nearly flat region of space because of Lorentz contraction along the beam- or z-axis. If 
the quarks and gluons stream freely into the surrounding space, then the distribution of 
particles would become locally highly anisotropic. In this case the stress tensor quickly 
acquires an oblate form with T^x ~ Tyy ^ T^z- The fastest growing mode of plasma 
instabilities then has its wavevector along the normal direction and generates a prolate 
contribution to the stress, which pushes the system towards greater isotropy. Theoretical 
discussions of thermalization typically start from the saturation picture of high-energy 
heavy- ion collisions, where there is initially a nonperturbatively large phase-space density 
~ of gluons with momentum of order the saturation scale Qs [H]. To be specific, 
we consider the extreme anisotropy case described by an effectively 5(j>z)-like initial 
distribution. 

^ For numerical simulations which take into account the backreaction of the soft fields on the hard 
particles using a Boltzmann- Vlasov treatment see [6]. 
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Fig. 1. LEFT: Fourier coefficients of the squared modulus of the gauge field versus time for three diff'erent 
momenta parallel to the z-axis. The curves correspond to the momenta in the legend from top to bottom. 
RIGHT; Primary and secondary growth rates for |A(t,p2:)p as a function of pz as measured on 96'^- and 
128^- lattices. (Taken from [S:.) 

The left graph of Fig. [1] shows the nonequihbrium time evolution of the color-averaged 
squared modulus of three different Fourier coefficients of the gauge field modes A{t^-p) for 
classical-statistical SU (2) gauge theory in three spatial dimensions. They are displayed 
as a function of time, normalized by the corresponding field values at initial timelfj All 
values are given in appropriate units of the initial energy density e. Here may be 

associated to particle numbers divided by frequency. The plotted low-momentum modes 
clearly show exponential growth starting at the very beginning of the simulation. In 
contrast to these "primary" instabilities operating at low momenta, one observes from 
the left graph of Fig. [T] that momentum modes at sufficiently high momenta do not 
grow initially. The higher wave number modes typically exhibit exponential growth at a 
"secondary" stage that sets in later, but with a significantly larger growth rate. The right 
of Fig. [1] displays the momentum dependence of the growth rates for obtained 
from a fit to an exponential, which is done separately for the primary and secondary 
growth rates. The secondaries arise from fluctuation effects induced by the growth in the 
lower momentum modes, which can be seen by taking into account (2PI) resummed loop 
diagrams beyond the hard-loop approximation 8 . This growth saturates when all loop 
diagrams become of order one, which leads to a highly non-linear subsequent evolution. 

In order to obtain an estimate in physical units one may consider an initial (Bjorken) 
energy density of about 5-25 GeV/fm'^ for RHIC experiments, and a projected factor of 
about two more for LHC energies. The inverse of the maximum primary growth rate for 
I Ait^'p) p is 7~ax — 4e, which then corresponds to a characteristic time scale of 

^ 1-2 - 1.8fm/c (RHIC) , 7"^ ^ 1-0 - 1.5fm/c (LHC). (4) 

One observes that the results are rather insensitive to the precise value of the initial 
energy density because of its scaling with the fourth root of e. For comparison, a time 
scale associated with the largest observed secondary growth rates is about a factor of 
three shorter than what is given in However, even though secondaries can reach 



^ The displayed results are for temporal-axial gauge. The characteristic growth rates have been verified 
also for Coulomb gauge [8]. 
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Fig. 2. LEFT: The ratio of modes of longitudinal pressure and transverse pressure as a function of time 
for different momenta [8]. RIGHT: Fourier coefficients of the squared modulus of the gauge field versus 
momentum for different snapshots of time. Fits to the spectrum at late time indicate a power-law with 
K = 3/2 for the classical-statistical theory. The latest shown time is t = 940e~^/'* (black squares) |17) . 

considerably higher growth rates than primaries, they start later. As a consequence, a 
certain range of higher momentum modes can 'catch up' with initially faster growing 
infrared modes before the exponential growth stops, as seen in Fig. [TJ This leads to a 
fast effective isotropization of a finite momentum range, while higher momentum modes 
do not isotropize on a time scale characterized by plasma instabilities. An example of a 
gauge invariant quantity, which reflects this behavior, is given by the Fourier modes of 
the stress tensor Tij. The left graph of Fig. [5] shows the ratio of modes for longitudinal 
pressure p) = T^^lt^p) with the corresponding transverse pressure for different 

absolute value of momentum. One observes a 'bottom-up isotropization' where for low- 
momentum modes the ratio turns to one rather quickly, while high-momentum modes 
only become isotropic at far later times. For the above energy densities one finds effective 
isotropization up to a characteristic momentum of about 



IpI <lGeV. 



(5) 



This is similar to the range of momenta for which hydrodynamic descriptions work, which 
typically show sizeable deviations from data above about 2 GeV. 

The subsequent slow evolution, after the fast exponential growth stopped, seems to 
provide a suitable condition for hydrodynamic models as effective low-energy descriptions 
based on sufficiently small gradients. What is the physics underlying this slowing-down? 
Most investigations concentrate on power cascades towards shorter wavelengths remi- 
niscent of Kolmogorov wave turbulence [21122117] . For the correlator of Fig. [T] such a 
behavior would be characterized by a power-law exponent k, with 



(6) 



The right graph of Fig. [2] shows the result of classical-statistical simulations in Coulomb 
gauge for transverse fields and momenta p aligned to the z-direction. One observes that 
as time proceeds more and more ultraviolet modes approach a power-law, which is rather 
accurately described by n = 3/2. For sufficiently high momenta this value may be com- 
pared to perturbative estimates. It has been pointed out that this value for k corresponds 
to the analytical result obtained from (2PI) resummed one- loop perturbation theory [17j . 



5 



Below we will see that the very same exponent is also found for inflaton dynamics. It 
has been noted that scaling behavior for gauge theories could deviate from scalar field 
dynamics and various estimates have been given for the value of k in the literature. In 
particular, the value two for k reported in Ref. [21] based on Vlasov equations seems ex- 
cluded by the x^~^^ for the classical-statistical simulation of Fig. [21 The difference might 
be attributed to the fact that the hard modes in the Vlasov treatment represent static 
sources that lead to a stationary-state solution, which is not present due to total energy 
conservation in the classical-statistical lattice gauge theory simulation. Simulations on 
much larger lattices might be required to see whether the value two for k could character- 
ize some intermediate-time behavior, which seems not the case in Fig. [2l Clearly, none of 
the classical approximations are sufficient to quantitatively address the late-time behav- 
ior, which should finally be characterized by a Bose-Einstein distribution for the gluons. 
While for low momenta the employed classical-statistical simulations are expected to give 
an accurate description for sufficiently high energy densities or occupation numbers, the 
high-momentum behavior will be altered by quantum corrections. It is remarkable that 
in the context of simpler scalar infiaton dynamics, where the quantum evolution can be 
addressed with present-day techniques, very similar phenomena are observed. 

3. Early universe reheating: A quantum example 

During inflation the universe exhibits a strongly accelerated expansion such that matter 
and radiation dilute very quickly. As a consequence, after inflation there is an enormous 
heating required in order to connect to the subsequent thermal history of the hot Fried- 
mann cosmology. In a large variety of infiaton models nonequilibrium instabilities play a 
crucial role for the process of thermalization. For instance, in chaotic inflationary models 
a parametric-resonance instability leads to an exponential growth of occupation num- 
bers |9|10|13j . In hybrid inflation, including D-term inflation, infiaton decay proceeds 
via a tachyonic (spinodal) instability of the inhomogeneous modes which accompany 
symmetry breaking [12ll4j . Despite the different underlying physical mechanisms, the 
subsequent evolution in these models after an instability follows very similar patterns. 

In quantum field theory dynamics of parametric resonance [13115] as well as tachy- 
onic preheating [T^ was studied in detail for A^-component infiaton models using (2PI) 
resummed 1/7V expansions to next-to-leading order (NLO) [23|24j . These type of approxi- 
mation are known to describe the late-time approach to thermal equilibrium characterized 
by Bose-Einstein [25123] or Fermi-Dirac distributions [2613] . respectively. Characteristic 
far-from-equilibrium phenomena, such as an early prethermalization of the equation of 
state, have been quantitatively studied in that context [3|27] . For a review see Ref. [28] . 

A characteristic example for chaotic infiation models is given in Fig. [3] which shows 
results for parametric-resonance amplification of occupation numbers in a weakly coupled 
0(4)-symmetric theory with quartic self-coupling A [13115] . The "occupation number" 

1 /2 

modes are defined as n{t,p) + 1/2 = ((!</?(<, p)p)c(|(/j(t, p)P)c) using connected cor- 
relation functions for the spatial Fourier modes (p{t, p) of the scalar infiaton field. The 
left graph shows the characteristic exponential early-time behavior of various momentum 
modes in units of the mass scale set by the rescaled initial infiaton field amplitude ~ 
{ip{t = 0,x.))y^X/6N. It is instructive to compare the initially fastest growing mode with 
momentum po and primary growth rate 270 with the behavior at a higher momentum 2po. 
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Fig. 3. LEFT: Exponential early-time behavior for the inflaton occupation number as a function of 
time for various momenta following a parametric-resonance instability il3i . RIGHT: Occupation number 
distribution for several later instances of time. At low momenta the nonperturbative critical behavior is 
rather well described by the infrared exponent k = 4 (dashed line) . For higher momenta the perturbative 
exponent k = 3/2 (dotted line) related to Kolmogorov wave turbulence is approached. (From 1151 .) 



The latter does not grow at first but exhibits exponential growth at a secondary stage 
with the larger growth rate 670, i.e. multiples of the primary growth rate. This occurs in 
the non-linear regime, where parametrically {\tp{t,p)\'^) ^ 0{N^^^\~^^^). Here occupied 
low-momentum modes act as sources for the secondary stage of enhanced amplification in 
a higher momentum range. The exponential growth stops when {\(p{t, p)P) ^ 0{N'^ ^~^) 
and the dynamics slows down considerably. At this nonperturbative stage all processes 
are of order unity and the system is strongly correlated despite the underlying weak 
coupling. It is striking to compare to the qualitatively similar primary and secondary 
growth stages for gauge fields in the context of plasma instabilities in Fig. [T] 

At later times the comparison between gluon and inflaton dynamics shows certain 
characteristic properties, which even quantitatively agree. The occupation number dis- 
tribution for the inflaton is displayed for several instances of time in the right graph 
of Fig. [3] [15] . The evolution is characterized by critical slowing down and the quasi- 
stationary solid curve exhibits different power-law regimes. At higher momenta the evo- 
lution approaches n(t, p) |p|~" with k = 3/2, which is indicated by the corresponding 
dotted line in Fig. [31 This is precisely the value characterizing the power-law behavior of 
gauge field fluctuations in the right graph of Fig. [21 The behavior is associated to weak 
Kolmogorov wave turbulence and can be obtained in perturbation theory. It should be 
emphasized that quantum corrections obstruct a power-law at very high momenta for 
which occupation numbers are of order one. 

The situation changes dramatically in the infrared below a momentum scale at which 
occupation numbers grow nonperturbatively large to n(t, p) ~ 247r/A for A ^ 1. A power- 
law behavior with strongly enhanced fluctuations emerges at low momenta according to 
the right graph of Fig. [31 This is indicated by the dashed line, which corresponds to a 
critical exponent k = 4. Analytically one finds that this is described by a nonthermal 
infrared fixed point solution of the renornialization group equations for the underly- 
ing quantum field theory, which is in the same universality class as the corresponding 
classical-statistical theory ^15.29 . It is still an open question whether gauge field theories 
exhibit this new type of infrared fixed points. The corresponding right graph of Fig. [2] for 
gauge field fluctuations in Coulomb gauge might indicate the slow emergence of a differ- 
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ent infrared regime [17j . and other indications seem to exist from reconsidering Vlasov 
results '21]. However, simulations for larger volumes and an ambitious analysis in terms 
of gauge invariant quantities such as stress-tensor correlation functions may be required 
to settle this important question. 

4. Conclusions 

Despite important differences the discussion of the thermalization process in heavy- 
ion collisions and cosmology after inflation shows remarkable similarities. In both cases 
nonequilibrium instabilities can lead to a fast period of exponential growth of occupation 
numbers. This is followed by a slow period, where the quantitative agreement concerns 
characteristic exponents or scaling functions describing scale-invariant properties of far- 
from-equilibrium dynamics. 

For a considerable class of inflaton models a detailed understanding in quantum field 
theory is available: Instabilities do not lead to fast thermalization, but they lead to an 
early prethermalization of bulk quantities such as the equation of state. The subsequent 
evolution approaches power-law behavior, which is described by nonthermal renormaliza- 
tion group fixed points of the underlying quantum theory. Similar to standard discussions 
in thermal equilibrium, one may distinguish between nonthermal infrared and ultraviolet 
fixed points. Important signatures of the nonperturbative infrared fixed points are related 
to strongly enhanced fluctuations on long-distance scales. This critical phenomenon is 
distinct from weak Kolmogorov wave turbulence relevant at shorter distances, which can 
be associated to the presence of an ultraviolet fixed point in the classical limit. 

Nonthermal infrared fixed points have the dramatic consequence that a diverging time 
scale exists far from equilibrium, which can prevent or substantially delay thermalization. 
They are approached from nonequilibrium instabilities without fine-tuning of parameters. 
Estimates of the reheating temperature after inflation have to take this nonperturbative 
physics into account. A conservative limit requires thermal equilibrium at a reheating 
temperature of order 10 MeV before Big Bang Nucleosynthesis. This can already rule out 
some simple weakly coupled inflaton models following parametric or tachyonic preheating 
dynamics, however, more realistic models have to be considered. 

In the context of heavy-ion collisions, QCD plasma instabilities may lead to a regime 
with scaling behavior. Classical simulations indicate possible fast "bottom- up" isotropiza- 
tion due to instabilities with characteristic time scale of 1-2 fm/c for low momenta of less 
than about a GeV. This appears to be near the range where hydrodynamic estimates 
give quantitative descriptions of experimental data. However, more reflned calculations 
in expanding geometries including quantum corrections for high-momentum modes may 
be required to settle this question. It is remarkable that simulations in classical-statistical 
lattice gauge theory show quantitatively the same scaling exponents as for the inflaton 
in the turbulent regime at higher momenta. The question of whether it is possible to flnd 
a nonperturbative infrared fixed point in QCD in possibly the same universality class is 
exciting. 

This work is supported in part by the BMBF grant 06DA267, and by the DFG under 
contract SFB634. Part of this work was done during the programme on "Nonequilib- 
rium Dynamics in Particle Physics and Cosmology" (2008) at the Kavli Institute for 
Theoretical Physics in Santa Barbara, supported by the NSF under grant PHY05-51164. 
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